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Abstract 

This paper presents a new approach to behavioral-social dynamics of 
human crowds. First order models are derived based on mass conserva¬ 
tion at the macroscopic scale, while methods of the kinetic theory are used 
to model the decisional process by which walking individuals select their 
velocity direction. Crowd heterogeneity is modeled by dividing the whole 
system into subsystems identified by different features. The passage from 
one subsystem to the other is induced by interactions. It is shown how 
heterogeneous individual behaviors can modify the collective dynamics, 
as well as how local unusual behaviors can propagate in the crowd. The 
paper also proposes a system approach to the modeling of the dynamics 
in complex venues, where individuals move through areas with different 
features. 

Keywords: Self-propelled particles, nonlinear interactions, crowd dynam¬ 
ics, mass conservation, active particles. 


1 Plan of the Paper 

The modeling of crowd dynamics can be developed, see m and the book 
at the three scales, namely microscopic (individual based), macroscopic (corre¬ 
sponding to the dynamics of mean averaged quantities), and to the intermediate 
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mesoscopic, corresponding to the dynamics of a probability distribution func¬ 
tion over the microscopic scale state of individuals. The latter approach is such 
that interactions are modeled at the micro-scale, while mean quantities, such as 
local number density and linear momentum, are obtained by velocity weighted 
moments of the aforesaid probability distribution. 

A critical analysis of the advantages and drawbacks of the different scales 
selected for the modeling approach are discussed in the review paper , where 
it is stated that the present state of the art does not yet allow well defined hall¬ 
marks to support an optimal choice. Therefore, due to these reasonings, a deep 
investigation has been developed to understand the complex dynamics at the 
microscopic scale, see, among others, [IH] and [23]. The result can contribute 
to implement both meso-scale models |4| or hybrid models [I], where the state 
of the system is defined in probability over the velocity direction and determin¬ 
istically over the velocity. A deep analysis of individual based interactions can 
contribute to derive, as we shall see, also hydrodynamic models mm- 

Macroscopic, hydrodynamic, models are of great interest in that they are 
far less computationally demanding than those at the other two scales. This 
requirement is particularly important when dealing with complex flows such as 
coupling pedestrian flows to vehicular traffic networks |13| . However, their main 
drawback is that the heterogeneous behavior of walkers gets lost in the averaging 
process needed to derive models, which ends up with hiding this important 
feature. An additional difficulty, well documented in the paper by Hughes ^5] . 
consists in modeling the process by which walkers select their velocity, namely 
direction and speed, in a crowd. 

In general, an important issue to be taken into account in the modeling 
approach is the nonlocal feature of interactions, as walkers are not classical 
particles and modify their velocity before encountering a wall. 

Finally, let us stress that challenging problem to be considered is the prop¬ 
agation of anomalous behaviors, which might be induced by panic conditions 
that can induce large deviations in the collective dynamics. Due to 
all aforementioned motivation the term social dynamics has been introduced in 
|7] to enhance the heterogeneous behavior of walkers, who might modify their 
strategy induced by interaction not only with other individuals, but also with 
the specific features of the environment where they walk. The interested reader 
can find in [21] a valuable investigation on hyperbolic equations generated by 
nonlocal interactions, see also [T4] . 

The present paper aims at tackling the aforesaid drawbacks in the case of 
first order models. These are simply obtained by mass conservation equations, 
which involve local density and mean velocity closed by phenomenological mod¬ 
els relating the local mean velocity to local density distribution. Although this 
approach is a simple way of looking at the dynamics, substantial developments 
are needed with respect to the existing literature, to achieve a realistic model¬ 
ing of the decision process by which individuals select the velocity direction and 
adjust their velocity to the local density conditions, as well as the heterogeneous 
behavior of the crowd. Moreover, this paper proposes a systems approach to 
model the dynamics in complex venues, where walkers move through ares with 
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different geometrical and qualitative features. 

The contents of the paper are presented in four sections. In detail, Section 
2 defines the mathematical structure underlying the modeling approach, which 
consists in a system of mass conservation equations for a crowd subdivided into 
various populations, which can be called after [^, functional subsystems. This 
structure acts as a background paradigm for the derivation of specific mod¬ 
els. Three structures are derived corresponding to a homogeneous population, 
an heterogeneous populations, but without social dynamics, and a population 
where social dynamics is accounted for. Section 3 shows how specific models 
can be derived according to the aforesaid structure. More specifically, three 
classes of models are proposed corresponding to the hierarchy of models defined 
in the preceding section. Section 4 presents some perspective ideas focusing on 
research plans, with more detail, the following topics are selected according to 
the authors’ bias: modeling the heterogeneous velocity distribution, propaga¬ 
tion of anomalous behaviors, as well as an introduction to a systems approach 
to model the dynamics in complex environments, where nonlocal interactions 
can play an important role in the passage through different areas. 


2 Mathematical Structures 

Let us consider the dynamics of a crowd in a domain E, which includes internal 
obstacles and inlet/outlet segments on the boundary 9S. This section searches 
for appropriate mathematical structures, which can provide the conceptual basis 
for the derivation of first order macroscopic models. It is assumed that the state 
of the system is described by local density and mean velocity to be viewed as 
dependent variables of the differential system. However, since we deal with first 
order models, it is also necessary looking for a functional relation linking the 
mean velocity to the density. 

Dimensionless quantities are used according to the following definitions: 

• p is the ratio between the number density n of individuals per unit area and 
the number um corresponding to the highest admissible packing density, 
namely p := 

• ^ is the dimensionless local mean velocity obtained by referring the di¬ 
mensional velocity v to the highest limit of the mean velocity VM^ which 
can be reached by a walker in a low density limit in optimal environmental 
conditions, namely ^ := 

Density and velocity depend on time and space coordinates, namely p = 
p(t,x.) =: p{t,x,y) and ^ = ^(t,x) =: ^{t,x,y). The independent variables also 
set by dimensionless quantities by dividing x and y by a characteristic dimension 
£ of S and time by Moreover, following [7], we introduce a number of 

dimensionless parameters that account for some specific features of the crowd 
already discussed in previous papers based on the kinetic theory approach: 
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• a S [0,1] models the quality of the environment where a = 1 stands for the 
optimal quality of the environment, which allows to reach high velocity, 
while a = 0 stands for the worst quality, which prevents the motion; 

• /3 S [0,1] models the attraction toward the direction of the highest density 
gradient, where /3 = 1 stands for highest attraction when all individuals 
follow what the other do, while /3 = 0 stands for the highest search of the 
less congested areas. 

These parameters are used in this section only at a formal level, while their 
physical meaning will be made precise in the next section devoted to derivation 
of specific models. 

Three levels of models are considered corresponding, respectively, to a crowd 
with homogeneous walking ability, a crowd subdivided into different populations 
featured by different walking abilities, and finally a social crowd where interac¬ 
tions can modify social behaviors, that can also include anomalous ones. 

• Homogeneous crowd: Let us first consider the derivation of the mass con¬ 
servation equation for an homogeneous crowd, where all individuals have the 
same walking ability. Classically, the said equation writes as follows: 

atP + Vx(p|) =0. (1) 

The closure of the equation can be obtained by modeling the dependence of 
4 on p by a phenomenological relation of the type ^ = ^[p](a,/3), so that the 
conservation equation formally writes as follows: 

5ip + Vx(p|[p](a,/3)) = 0, (2) 

where square brackets denote that functional, rather than functions, relations 
can be used to link the local mean velocity to the local density. For instance, 
walkers select their preferred directions based not only on the local density, but 
also on density gradients. Therefore, specific models can be obtained by some 
heuristic interpretations of physical reality leading to ^ = ^[p] and inserting 
such model into Eq.([2]), which is the mathematical structure to be used for the 
modeling. 

• Heterogeneous crowd: Let us now consider a more general framework, 
where walkers are subdivided into a number n of populations, the labeled by 
the subscript i, corresponding to different levels of expressing their walking 
abilities. Therefore, the state of the system is defined by a set of dimensionless 
number densities 


n 

Pi=p^{t,x.), i = p{t,x) ='^pi{t,x), (3) 

i=l 

The subscripts correspond to a discrete variable, modeling the walking ability, 
with values corresponding to the the subscripts i = 1,... ,n, being i = 1 and 
i = n, respectively, to the lowest and highest motility. 
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The new structure simply needs the following modification: 

+ Vx (p*|j[p](a,/3)) = 0, = (4) 

where the modeling of the mean velocity differs for each population [p]. 

Therefore, the structure consists in a system on n nonlinear PDEs equations. 

• Social crowd: Let us finally consider a dynamics, where each individual 
expresses a certain social variable, for instance (but not only) panic, with values 
in the interval [ 0 , 1 ], where the extremes denote the lowest and highest level of 
social expression corresponding to such variable. Additional notations and a 
substantial development of the structure Eq.(|3]) are needed. 

1. A social variable, say w, is introduced with discrete values j = 1,.. .,m, 
where j = 1 denotes the lowest level of expression, e.g. absence of panic 
w = 0, while j = m denotes the highest admissible value w = 1 . 

2. The overall system is subdivided into N = n x m populations labeled by 
the subscript ij. 

3. The density in each population is denoted by pij(t,x), therefore the total 
density is 

n 771 

i=i i=i 

4. A source term Sij[p] models the transition from one functional subsystem 
to the other due to social interactions. The formal expression of such term 
is as follows: 


Szj = Sij[p]{a,l3,-i), with EE 5'y[p](Q!,/3,7) = 0, (5) 

i=i i=i 

where 7 is a parameter triggering the transition into the ij-state, which 
depends on the specific type of social exchange. 

The mathematical structure now consists in a system oi n x m nonlinear 
PDEs as follows: 


^tp^J + \'^{p^J^ij[p]ia,P)) = S^j[p]{a,l3,-f) , i=l,...,n j = l,...,m, ( 6 ) 

and where the modeling of the mean velocities = ^ij[p] differ for each ij- 
population. 

It is useful, for sake of completeness, considering also the case of a crowd 
with a homogeneous mobility, but heterogeneous social state. Such a structure 
can be useful to study the propagation of anomalous behaviors. A simplification 
of Eq. (| 6 ]) yields: 

dtPj+ V^{pj£,j[p]{a,l3)) = Sj[p]{a,P,-i), j = l,...,m. (7) 
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Remark 2.1 Structures defined by Eqs. and (2.6) constitute the under¬ 

lying framework to derive specific models. In all cases, the derivation depends 
on the environment, where the crowd moves, not only its quality, but also its 
shape. In fact, the model requires implementing nonlocal interaction with the 
walls. Moreover, transferring the structure W into a model requires a detailed 
analysis of the social dynamics. 

Remark 2.2 Specific models, where some social dynamics is taken into account, 
need teehnieally a specification of the kind of social information exchanged by 
walkers. 

Remark 2.3 The formal structures proposed in this section need three types of 
parameters, namely a, j5, and 7 , eorresponding to three specifie features of the 
dynamics. Their practical implementation into models will be clarified in the 
next section. 


3 Derivation of Models 

Let us now consider the derivation of specific models consistent with the frame¬ 
works presented in the preceding section. Therefore, three different classes of 
models are derived in the following subsections according to conceivable require¬ 
ments of applications. More in detail, we consider both homogeneous and het¬ 
erogeneous populations in absence of social interactions, and an heterogeneous 
population when social exchanges take place so that individuals are allowed to 
shift from one subpopulation to the other. The study of panic conditions or, 
more in general, propagation of anomalous behaviors can be treated in the latter 
case. These contents are presented in the next subsection, while the last one 
proposes a critical analysis looking ahead to research perspectives. 

3.1 Modeling an homogeneous crowds 

Let us consider the simple case of an homogeneous crowd, where all individuals 
behave in the same manner. Derivation of models requires simply to describe 
analytically the dependence of ^ on the local density distribution. Polar coor¬ 
dinates are used for the velocity variable, so that 

4 ^ (cos0i-I-sin0j) := (8) 

where 6 is the angle between ^ and the a;-axis of an orthogonal system in the 
plane, i and j are the unit vector of the aforesaid axes, ( is the velocity modulus 
of ^ occasionally called speed. 

The idea that pedestrians adjust their dynamics according to a decision 
process, which can modeled by theoretical tools of game theory |31] . was already 
introduced in [3] and subsequently developed in a sequel of papers mil [7]. 

In detail, the modeling proposed in this paper takes advantage of the ap¬ 
proach to individual based interaction proposed in [7] , where the decision process 
is supposed to act according to the following sequence: 
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1. Walkers move along the direction us, forming an angle 9 with respect to 
the x-axis, selected according to a decision process which account the 
following trends: search of the exit, avoid walls, search of less congested 
areas, and attraction toward the density gradients. Details are given in 
the following. 

2. Once the direction uj has been selected, pedestrians adjust their velocity 
modulus to the so-called local perceived density pg along lj. This quantity 
is defined, according to [5] as follows: 

Pe = Pe[p]=P+ .y [{^ - p) H(dep) + pH{-dgp)] , (9) 

^l + [dgpY 

where dg denotes the derivative along the direction 0, while H{-) is the 
Heaviside function i?(- > 0) = 1, and iJ(- < 0) = 0. Therefore, posi¬ 
tive gradients increase the perceived density up to the limit p = 1, while 
negative gradients decrease it down to the limit p = 0 in a way that 

dgp oo => Pe 1 , dgp = 0^ pI = p, dgp -oo => Pe -)■ 0 . 

( 10 ) 

3. The relation ^ = ^(pe; a) between ^ and Pe depends on the quality of the 
environment. In the attempt of reproducing empirical data [m [Ml [501131 
[33], is modeled by a polynomial of pg fulfilling the following constraints: 
^(0) = a; ^'(0) = 0; ^(1) = a; ^'(1) = 0, where prime denotes derivative 
with respect to pg. 


The first step of the above process leads, according to [7] to the following 
expression of the preferred direction: 


f ^ (1 - p)l^T + p [Pl^S + (1 - Id^v] 


( 11 ) 


where 

while the parameter /3 may account for panic conditions [7] . Hence, cos 9 = 
U3 ■ i,sin 6 * = a; • j. 

Neglecting the influence of the walls, the second and third step yield the 
following polynomial expression: 


Cf = ■Cf’M = a(l - 3pf-f 2p;^)[p]. (13) 


These two quantities can be inserted into the framework Eq. o, which 
writes as follows: 

9tP + Vx(ap(l-3pf+ 2pf)u;f(p,x;)) =0. (14) 
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The role of the nonlocal interactions with the wall and obstacles can be 
modeled by adapting the hallmarks of [7] to this present case. The strategy 
pursued is that the aforementioned interactions, which replace the classical local 
boundary conditions, define a preferred walking direction according to a two- 
steps procedure. As a first step, walkers change the direction of motion to ujp. 
As a second step, walkers may further change the direction from cjp to u: if two 
conditions occur: 

1. The walker’s distance from the wall, d, is within a given cutoff distance, 

dw 5 

2 . the walker’s velocity is directed toward a wall. 

In detail, if ^ <0 and d < d^ , the components oi ujp are decomposed 

into the normal and tangential components, respectively 

a;i(p,x;/3) = (n(g>n)a;F(p,x;/3), and W2(/9, x;/3) = (I - n 0 n)wF(p, x;/3), 

(15) 

where n is the unit vector orthogonal to the wall. Subsequently, it is assumed 
that the velocity component normal to the wall decreases linearly while ap¬ 
proaching the wall and becomes naught in the limiting case of a pedestrian in 
contact with the wall, hence 

d ( d? \ 

‘^(p,x;/3) = —a;i(/9,x;/3)n-H ( 1 - —wJ(p,x;/3)J r, (16) 

where r is the unit vector tangential to the wall. The model is obtained by 
Eq. (fT^ simply by replacing with us. 

3.2 Heterogeneous crowd in absence of social dynamics 

The reference mathematical structure is now given by Eq. o, which consists 
in a system of n PDEs. Hence, the modeling problem consists in modeling 
the mean velocity for each population of pedestrians, which, as an example, 
can correspond to slow, mean, and fast individuals. Hereinafter, the model 
is presented in the case of a system in absence of boundaries, referring the 
technical generalization to the approach proposed in the preceding subsection. 
By assuming that pedestrians exploit the quality of the environment depending 
also on the quality of their walking ability, the following model is obtained: 

^^ = ap\ = o^'-{^-^pf+ ‘^pf)[p] (17) 

which inserted into the structure yields: 

+ Vx ^ (1 - -b 2pg^)[p]a;j(p.x;/3)^ = 0, i = l,...,n, (18) 

where the preferred direction is now computed for all components of the fluid 
viewed as a mixture. 
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3.3 Modeling crowds with social dynamics 

Let us now consider the dynamics of individuals according to the general frame¬ 
work, which introduces a discrete social variable. Specific models can be ob¬ 
tained by a detailed description of the source terms 5^. Therefore, let us 
consider the modeling of transitions, induced by interactions, across functional 
subsystems to the other. These can be viewed as net flows, namely the inflows 
into the j-subsystem from the {i — 1 )- and {i + l)-subsystems, and minus the 
outflow from the j-subsystem into the {i — 1) and {i + 1) subsystems. A formal 
expression is as follows: 

S^3 [P] = Sfj [p] - S-j [p] = (7+ Pi(j-i)Pij - l~p^(:j+i)Pij) 

~ {li PiU+l)Pij ~ li (19) 

where and 7 ” are positive constants modeling, for each i, the flow from 
the lower to the upper state and vice versa and where the flow depends on 
the density of the state origin of the flow and that attracting it. Letting now 
7 i = 7 / — 7 ~, which is a parameter with positive or negative values, yields 

S^3[p\ = Stj[p\ - S-j[p\ = 7 * iPi(j-i) - PiU+i))- (20) 


Therefore 

e > 0 : Pi -1 < pi+i ^ Sij < 0 , Pi -1 > pi+i => Sij > 0 . ( 21 ) 

Remark 3.1 This result can he specifically related to modeling panic, which 
is a situation such that initially Pi(j-i) — Pij, while increasing of panic shifts 
individuals to high values off, namely Pi{j-i) < Pij- This topic will be further 
discussed in the last section. 


4 From a Critical Analysis to Perspectives 

Three classes of models have been proposed in preceding sections, in sequence, 
corresponding to: A population with homogeneous behavior, an heterogeneous 
walking ability and strategy, and to a dynamics which involves social behav¬ 
iors, specifically panic. The advantage of the approach proposed in this paper 
mainly relies on the simplicity of the mathematical structures used for the mod¬ 
eling mass conservation closed by heuristic models linking mean velocity to local 
density distribution. This strategy leads to models which present a lower com¬ 
putational complexity with respect to the kinetic type approach [7]. On the 
other hand some descriptive ability is lost in the averaging process induced by 
the hydrodynamic approach. 

One of the aims of this paper has been developing a new approach suitable 
to include, as far as possible, some features of the heterogeneity, which appears 
in crowd dynamics. This final section aims at understanding how far the math¬ 
ematical structures proposed in the preceding sections can be further developed 
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to model additional complexity features, heterogeneity in particular, already 
described in [7]. This development takes advantage of Gromov’s hint toward a 
deep investigation of mathematical structures with the aim of discovering their 
richness [53]. 

The contents focus on the following topics selected, among various ones, 
according to the authors bias: Additional reasonings on modeling different types 
of social behaviors, modeling heterogeneity in the velocity distribution, and 
finally a system approach to the dynamics in complex venues, where the passage 
to different types of areas needs to be carefully modeled by taking into account 
nonlocal interactions. These topics are treated in the next subsections, which 
although not formalized by mathematical equations, propose various hints to 
tackle the related conceptual difficulties. 

4.1 Modeling different types of social behaviors 

The modeling proposed in Subsection 3.3 is focused on panic conditions. How¬ 
ever, a variety of different types of exchanges can appear. The modeling of 
interactions can be achieved by suitable development of theoretical tools of 
game theory Eiisi]. A recent literature on the mathematics of social systems 
provides useful hints toward this approach, see among others HlHKinillD]. In 
particular, the following types of social exchanges can be considered: 

• Consensus dynamics: When walkers approach, during interactions, their 
mutual social status, namely individuals with low values of the activity 
variable shift to the higher value, while those with an high value shift to 
the lower one. 

• Competitive dynamics: When individuals behave in the opposite way and 
increase the distance between their mutual social status, namely individ¬ 
uals with low values of the activity variable shift to a lower value, while 
those with an high value shift to an higher one. 

• Leaders dynamics: When the presence of leaders induces a different dy¬ 
namics as leaders do not change their status, while the other individuals 
shift toward the status of the leaders. 

A detailed formalization of this qualitative description leads to the modeling 
of the fluxes, which appear in Eq. (2.6) and hence to mathematical models 
related to the three different types of dynamics under consideration. 

4.2 Modeling heterogeneity in the velocity distribution 

The modeling approach proposed in this paper shows a drawback typical of all 
hydrodynamic models, namely the averaging process hides the heterogeneity 
features of the system. In this present paper, heterogeneity of individuals is 
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taken into account by the subdivision into functional subsystem, but the het¬ 
erogeneity in the velocity distribution is lost. This feature obliges one to use 
empirical data to refer the speed to the local density distribution. 

The empirical relation defined in Eq. m links the speed to the parameter a 
and to a polynomial expression of the perceived density, where a has the mean¬ 
ing of highest observed speed within a representation which is one-dimensional 
due to the fact the velocity direction has been selected according to the decision 
process followed by walkers. A more refined approach might consider ^ as the 
first order moment of a probability distribution, evolving in time according to 
interactions at the micro-scale to be properly investigated. 

This objective can be achieved by taking advantage of the modeling approach 
proposed in [12] for vehicular traffic according to the following hallmarks: 

• Walkers select their optimal velocity direction oj according to the deter¬ 
ministic decision process proposed in the preceding section; 

• Speed is assumed to be given by a random variable u(t; a, x, p) depending 
locally on a and on the perceived density; 

• The variable v is linked to a probability distribution P(t, v\a, x, p), whose 
dynamics is derived, as in [12] , by individual based interactions depending 
on a and on the perceived density; 

• The mean speed ^ is obtained as first moment of P{t,v) for a system, 
where the overall dynamics is delivered by coupling the dynamics of p to 
that of P{t, v). 
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